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ABSTRACT 

Recent data analysis of polar motion indicates the presence of a component 
with periodicity corresponding to the motion of the lunar ascending node. An 
investigation of the tidal response of the ocean to long period forcing functions 
has been conducted. The results of the investigation indicate the possibility of 
excitation of a wobble component with the amplitude and frequency indicated by 
the data. An enhancement function for the equilibrium tide has been postulated 
in the form of an expansion in zonal harmonics and the coefficients of such an 
expansion have been estimated so as to obtain polar motion components of the 
required magnitude. 
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THE ENHANCED NODAL EQUILIBRIUM OCEAN 
TIDE AND POLAR MOTION 


I. INTRODUCTION 

Recent data analysis by Markowitz (1979) indicate the existence of a wobble component with 
periodicity and phase corresponding to the motion of the lunar ascending node. 

The objecdve of this investigation is to ascertain the possibility of the existence of a 18.6 
year wobble component due to a modified equilibrium ocean tide. The study of the equilibrium 
response of the oceans dates back to Darwin (1886). More recently Proudman (1960) reached the 
conclusion that the tidal constituent with nodal period will follow the equilibrium law while the 
semiannual and annual constituents will probably follow it. 

Wunsch (1967) tested the equilibrium hypothesis by computing periodograms for the fort- 
nightly and monthly tides on islands of the Pacific which he found to deviate significantly from 
equilibrium. He considers the nodal tide to tend towards equilibrium with a certain degree of 
admissible uncertainty. 

The ocean responds not only to the gravitational potential of the Moon and the Sun, but also 
to the second degree potential of the Earth rotation. The “pole tide" or ocean response to the 
Chandler wobble has been analyzed by various investigators. Haubrich and Munk (1959) analyzed 
mean monthly values of sea level from 1 1 tide stations and found the average pole tide with period 
of 14 months to have an amplitude twice that predicted by equilibrium theory. 

Hosoyama, ct at., (1976) found the amplitude ratios of the observed to equilibrium pole 
tides to increase at higli and low latitudes in the northern hemisphere; they also find latitude 
dependent phase delays and advances with implications concerning the possible excitement of the 
Cliandler wobble. Naito (1977, 1979) computed the secular variations of the amplitude and phase 
of the observed pole tide for the period 1900-1964 and compared them with those of the 
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equilibrium pole tide, he found the largest amplitudes of the observed pole tide to take place in 
the coasts of the Baltic and North Seas, he concludes that the observed pole tides seem to have a 
certain relation with the equilibrium tides, but to have their own secular variations. Dickman 
(1979) studied the effects of the pole tide on the (Thandler wobble assuming pole tide amplitudes 
up to 10 times equilibrium both for the case of global enhancement and for regional enhance* 
ments in the North and Baltic seas. Dahlen (1976) developed a general theory to determine the 
influence of the pole tide upon rotation of the Earth. 

The concept of an enhancement function which modifies the equilibrium tide can be applied 
also to the ocean response to the gravitational disturbing forces. Such an enhancement function 
can be expressed in terms of spherical harmonics with adjustable coefficients. It is then possible 
to express the tidal contributions to the products of inertia in terms of the enhancement function 
coefficients which can be estimated so as to obtain the polar motion components indicated by the 
analysis of the data. The estimated coefficients then can be used to predict the modified behavior 
of the equilibrium tide. 

2. SOLUTION TO THE LIOUVILLE EQUATIONS 

Tlie Liouville equations of motion were first given by Liouville (1858). The following 
assumptions are now made, 

(1) tlie external moments and relative angular momentum terms vanish. 

(2) the moments of inertia are constant and considerably larger than the products of 
inertia; the equatorial moments of inertia are equal. 

(3) the component of angular velocity is a constant and much larger than co^ and cOy. 

Neglecting products of small quantities the equations of motion then become 

( 2 . 1 ) 

-iy^CJ^ + AWy - 1x2^2 + (A - C) = 0 
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Now let the products of Inertii be given by 


*xt " *xs ♦ *xt 
V* " *ys * ^yt 


(2.2) 


where the superscript r denotes the contribution due to rotitional deformstion of the solid Earth 
and T indicates the contribution due to the ocean tide. The rotational deformation is known to 
be given by 


r * 

>XI “ - J kj — WxW* 

*yi " "7‘'2 ^ 


(2.3) 


where R denotes the radius of the Earth, G is the gravitational constant and k} is the second degree 
Love number. 


Equations (2.1) can then be written as follows, 
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where 


Wx + 


(A + a,) 
1 


ih “ nfy) 


^ a • 

civ + n^cjv * ■ (fv + nf-) 

^ ^ (A + a,) ^ 


I RS , 


KC - A) - 1 

n « to, 

(A + a,) 

(x “ ” ^yz^i 


(2.4) 


(2.5) 
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lyx " ^yz “ ^yx^ 

Then 

Wx “ cos({t - 



Wy * Wy C08(ft - ^<jy) 


M 

II 

{k? + + 2 K,K 2 Sin(«xz - 6 yz>] 

^t/2 

Wy = 

+ KJ + 2K3K4sin(^yx - 0xz)] 


* arc tan 

•^K,sin0xz KjCos^yxVfKiCOs^xz 

- Kjsin^yj) 

= arc tan 

^K^sin^y^ ^ R4 ^^s0xzV(K 3 ^^^^y z 

- Ka'in^xz) 

Ki = 

-MxzWz(f* + nwjI/KA + a,Kn* - 

f*)l 

II 

MyzWzf(**^z n)/UA + a,Xn^ - 

f*)l 


K, = (Myj/Mxz)K, 

K4 » -(Mxz/My^)K2 


The polar motion components are given by 

X = (<4)x/(^z)K 
y * (Wy/u»j)R 


( 2 . 6 ) 


(2.7) 


( 2 . 8 ) 


3. EQUILIBRIUM OCEAN TIDE 

The subject of an equilibrium ocean tide including the effects of ocean loading and the self 
attraction of the water has been the object of various investigations; Hcndersholt (1972), Dahlen 
(1976), Agnew and Farrell (19'^8) and others. It follows that a global equilibrium tide is given by 
the following expression: 


£ = 


1 - ®n(^ ® 

3 P 




2n + 1 p,; 


(3.1) 
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wher« U„ is the disturbing potential due to the mass of the tide generating body and its motion 
relative to the Earth; k„, h„, k|, and h|| are the Love numben of degree n, p denotes the density 
of water and p^ is the mean density of the solid Earth, g is the acceleration of gravity. 


in order to introduce the concept of a modified equilibrium tide an enhancement (Unction is 
postulated and its functional form is assumed to be given by an expansion in spherical harmonics 
of the following form, 

E(0) = r§P§ + + ^5*1 (3.2) 

where the P’s arc the Legendre polynomials and the r’s are undetermined coefficients. The modified 
equilibrium tide is then given by 

(m = (3.3) 

Tlie uninodificd tide is recovered when 

rg = 1. r? = r§ = rg = 0. 


Now let SI denote the surface area covered by the oceans and define as foUows: 

{n(ds (3.4) 

where the surface integral is taken over the area of the continents. represents the quantity 
which must be added to in order to satisfy conservation of mass. Note that. 


ff ~ ff ■ ff 

''•'C •'•'Sohere •^•'Oceans 

fL 


^Sphere 

“ ff 

^Oceans JJSphtte 


Oceans 
f(0.^)tM dS 


(3.5) 


= PS’ 


n m 


a^cosm\Jf’ 
. b{” sinm^ 


■{; 


over continents 
over oceans 


(3.6) 


consider the ease when the disturbing potential is given by 


(3.7) 

qj being a function of the position and mass of the disturbing body. Making use of Equations 
(3.1)-(3.7): 



h 


(1 + k, - h,) 

1 - + k'j - h',) 



g 


(3.8) 


The complete expression for the modified equilibrium tide is then 


(3.9) 


4. PRODUCTS OF INFRTIA 

I'he contributions to the products of inertia due to a modified equilibrium tide are given by 


Sptierc 


■ 7 / 

Sj. 

■ 7 / 

J J Sphere 


f(d, 1 ^) xz dm 


dm 


(4.1) 


where 

.'f - RsindcosJ/ 
y == Ksintfsin^ 
z = R cosO 


() 


dm - pIj;,(9)(Rdd • Riindd^] 


Making use of Equations (3.1H3<9) yields 



Ll75\b}/ \1225 1575/ \bj/ \24255 99099/ \b|/ 

165165 \b^/J 5 Vbi/Bjj 


(4.2) 


The integration of tlie product of three and four spherical harmonics has been performed by means 
of the 3-j symbols of Wigner (Rotenberg, et al., 1959; Winch and James, 1973). 

5. NUMERICAL RESULTS 

The expressions for the products of inertia given by Equations (4.2) can be put in the form 
of Equations (2.6) since 

q§ = “(Gd/R 2) 53 MiCOs(o,) (5.1) 

i 

where Gp is Doodson’s constant and M{, Oj are the amplitudes and arguments for the various tidal 
constituents. The principal terms of the low frequency tides ai'e given by Cartwright and Edden 
(1973). For the purpose of this investigation the nodal term is the one of interest, then 

M » -0.06556 

a = ft - 259.18328° (5.2) 

f = 2ir/I8.613 years 
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q>odi. i89l J^ecember 31, 12h, Om, Ot 
ephemetb time 

The tolution of the problem con^su in eetimeUnt the vtluee of r^,rf ,rj| endf^ 4>peuint in 
Equation (4^) which will yield polar motion components aatiifyinf the remits from data analysis. 

Msriiowitz (1979) analysed 79 years of ILS (International LaUtude Service) data and 17 
years of IPMS (International Polar Motion Service) data, he obtained the foUowing remits: 

ILS 

X “ (28 ± 13) cos (ft - 1903.7 ± 1.2y) 

y » (22 ± 13) cos (ft - 1904.4 ± l.5y) 

IPMS 

X “ (22± 13)cos(ft - 1905.7 ± 1.5y) 

y = (25 ± 13) cos (ft - 1906.7 ± 1.4y) 

Tlie results are expressed in a geodetic coordinate system in cm units. The polar motion 
components corresponding to the case of a non’^odified equilibrium tide are obtained by letting 
r§ = 1 , r? = r§ « r^ * 0, this yields 

X » 0.31 cos (ft - 1906,14) 
y = 3.02 cos (ft - 1913.38) 

In order to estimate the values of the coefficients that will tit the polar motion (lata a general 
purpose adaptive iterator for nonlinear problems (Campbell, rl al.) has been used. The results are 
given below: 

lt§ 

rg » ' .1260094 

r\’ » 2.7151830 
rg » 3.1766850 
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r5 - -0.92669194 
X • 27.48 cos (ft - 1905.37) 
y « 22.84 cos (ft - 1904.97) 

IPMS 

rg - -2.0466899 
r? • 1.7925059 
> 2.0108450 
rg » 0.0091169942 
X ■ 2 1.27 cos (ft - 1905.44) 
y « 25.85 cos (ft - 1905.05) 

In both cases the cocffldents were estimated so as to obtain polar motion components within 
1 cm of tlie mean values given by Markowitz analysis. No constraints were imposed on the values 
obtained for the phase angles. 

The corresponding enhancement functions and tide heights follow from Equations (3.2) and 
(3.9), they are shown in Figures (5.i)-(S.5) below, also shown is the unmodified equilibrium tide. 

6. CONCLUSIONS 

The results indicate that a modifled equilibrium tide could provide the excitation required to 
generate the polar motion component detected by Markowitz. However since no ocean dynamics 
have been incorporated into the formulation of the enhancement function the results have to remain 
speculative. The concept of enhancement is not a new one althou)^ its application has been limited 
to the pole tide, in such cases different investir'^tors have considered the possibility of enhancement 
reaching values up to 10 times the equilibrium, also analysis of tidal data has indicated the existence 
of latitude dependencies as well as magnifications in shallow seas. In that light and by c'*<nparison the 
nuignitude of enhancement obtained in this study does not appear altogether exorbitant, especially in 
the case of the IPMS data. Neverthelesr the results should be considered as indicative of a possibility 
rather than as a quantitative determination of ocean behavior. 
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APPENDIX 1 


Numerical values of the constants used in the investigation. 


Moments of inertia. 


A “ 8.016(10^^)gm - cm^ 

C » 8.043(10^^ )gm - cm* 
Wj “ 2 ir/86400rad/8ec 
R = 6.378(10«)cm 
kj = 0.30 

G = 6.67(10'®)cm* dynes/gm* 


With the values above and making use of Equation (2.S) one obtains a 44S day Qiandier period. 


(1 + kj - hi) 

1 - Old + k^ - hi) 


= 0.86465 


p = 1.03gm/cm* 
g = 980 cm/sec* 


The values of the coefflcients a||* and bj|' appearing in the ocean function, Equation (3.6) 
are those given by Balmino, Lambeck and Kaula (1973). The proper normalization factors have 
been applied in order to maintain consistency throughout the calculations. 


Gp = 2.627723(10^)cm*/sec* 




FIGURE CAPTIONS 


Figure 5.1 Unmodified Tide Height 
Figure S. 2 Modified Tide Height (ILS) 
Figure S.3 Modified Tide Height (IPMS) 
Figure 5.4 Enhancement Function (ILS) 
Figure 5.5 Enhancement Function (IPMS) 






